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NONLINEAR  THEORY  AND  EXPERIMENTAL  OBSERVATIONS 
OF  THE  LOCAL  COLLISIONAL  RAYLEIGH-TAYLOR  INSTABILITY 
IN  A  DESCENDING  EQUATORIAL  SPREAD-F  IONOSPHERE 


INTRODUCTION 

Recently,  coincident  rocket  and  radar  observations  of  equatorial  spread- 
F  were  made  near  Kwajalein  in  the  Marshall  Islands  [Szuszczewicz  et  al. ,  1980 
Tsunoda ,  1980a].  The  in  situ  rocket  data  showed  major  plasma  depletions  dis¬ 
tributed  throughout  a  downward  moving  F-layer  with  regions  of  smaller  scale 
in  situ  irregularities  tending  to  be  collocated  with  positive  gradients  of 
electron  density.  The  VHF  radar  measurements  showed  that  the  backscatter 
"plumes",  in  addition  to  having  large  vertical  extent  from  the  bottomside  to 
the  topside  [see  also  Kelley  et  al.,  1976;  McClure  et  al. ,  1977;  Woodman  and 
LaHoz ,  1976J,  were  identifiable  with  the  plasma  depletions  [Tsunoda  and  Towle 
1979;  Tsunoda ,  1980b]  and  could  be  characterized  by  a  growth  and  decay  phase 
L Tsunoda,  1980aJ.  The  recent  rocket/radar  comparison  [Szuszczewicz  et  al., 
1980]  specifically  showed  that  the  most  intense  radar  returns  in  the  topside 
F-region  were  collocated  with  the  upper  region  of  a  decay-phase  depletion. 

Several  features  of  the  small  scale  irregularity  3m)  signatures 

seen  by  the  radar  backscatter  can  be  explained  by  various  plasma  kinetic 

instabilities  [for  a  review,  see  Ossakow,  1979]  which  are  presumably  driven 

by  steep  gradients  within  developing  plasma  bubbles.  On  the  other  hand, 

much  evidence  now  exists  that  the  large  scale  irregularities  result  from  a 

Rayleigh-Taylor  instability  [Balsley  et  al . ,  1972;  Haerendel,  1974]  driven 

by  the  ambient  bottomside  plasma  density  gradient.  Global  [Scannapieco  and 

Ossakow,  1976;  Ossakow  et  al. ,  1979;  Zalesak  and  Ossakow,  1980]  and  local 

[Keskinen  et  al..  1980]  numerical  simulations  of  the  collisional  Rayleigh- 

Taylor  instability  have  successfully  reproduced  the  large  plasma  density 

depletions  (bubbles)  and  associated  spatial  nower  spectra  that  have  been 

experimentally  observed. 

Manuscript  submitted  September  30,  1980. 


The  most  recent  in  situ  rocket  observations  [ Szuszczewicz  et  al.,  1980] 


were  made  in  an  F-layer  that  was  moving  downward,  implying  the  existence  of 
a  large  scale  westward  electric  field.  The  importance  of  this  electric  field 
in  describing  the  development  of  plasma  bubbles  has  been  discussed  by  several 
authors  [Ossakow  and  Chaturvedl,  1978;  Ossakow  et  al.,  1979;  Ott ,  1978; 
Anderson  and  Haerendel,  1979].  However,  all  previous  analytical  and  numerical 
studies  of  the  nonlinear  collisional  Rayleigh- Taylor  instability  have  not 
included  the  effects  of  an  upward  or  downward  moving  F-layer  in  realistic 
geometries,  nor  have  they  had  the  opportunity  for  an  almost  ideal  comparison 
with  experimental  results. 

In  this  paper  we  study  the  nonlinear  evolution  of  the  collisional 
Rayleigh-Taylor  instability  in  a  vertically  convecting  F  layer  using  both 
analytical  and  numerical  techniques.  In  addition  we  compare  the  simulations 
with  the  recent  in  situ  rocket  data  which  provided  direct  measurements  of 
gradient  scale  lengths  and  density  fluctuation  power  spectra. 

EXPERIMENTAL  RESULTS 

By  2100  (LT)  on  17  July  1979  the  bottomside  of  the  F-region  at  Kwajalein 
had  risen  to  an  approximate  altitude  of  400  km.  The  F-layer  then  began  a 
downward  drift  with  an  approximate  velocity  VD  =-  10  m/sec  [Szuszcezewicz 
>  198°;  R-  Tsunoda,  private  communication,  1980]  and  a  simultaneous 
onset  of  spread-F.  With  the  F-layer  drifting  downward  and  spread-F  con¬ 
ditions  continuing,  a  rocket  (designated  PLUMEX  I)  was  launched  (00:31:30, 

July  17,  1979  LT)  when  the  bottomside  F-layer  had  descended  to  an  altitude 
below  300  km. 
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A  pair  of  pulsed  probes  [Szuszczewicz  and  Ho lines,  1977;  and  associated 


references]  were  diametrically  extended  from  the  forward  end  of  the  rocket 

payload.  The  sensing  elements,  constructed  from  tungsten  wire,  were  isolated 

from  their  extension  booms  by  coaxial  guard  electrodes  driven  at  the  same 

potential  as  the  probes  themselves.  One  of  the  probes,  defined  as  the 

I-probe,  operated  with  a  baseline  voltage  VR  -  lv,  yielding  net  ion  baseline 

current  1^.  The  other  probe,  defined  as  the  E-probe,  operated  with  V  +  2v, 

D  B 

0 

yielding  net  electron  baseline  current  I  .  Both  probes  generated  complete 

b 

current-voltage  characteristics  in  T  «=  400  msec,  yielding  absolute  values  of 

s 

and  Tg  at  an  approximate  2.5  Hz  rate.  Maximum  I  sampling  occurred  at 

2048  Hz,  resulting  in  0.5  meter  spatial  resolution  for  relative  electron 

density  fluctuations  at  a  vehicle  velocity  of  1  km/sec.  (See  Szuszczewicz 

and  Holmes,  1980  for  instrument  details  and  complete  payload  configuration.) 

Figure  1  displays  the  upleg  measurements  of  relative  and  absolute 

electron  density  as  extracted  from  the  pulsed-plasma-probe  data.  Absolute 

values  of  electron  density  were  determined  by  conventional  analyses  of 

Langmuir  probe  characteristics  [Chen,  1965]  with  appropriate  care  to  eliminate 

perturbing  effects  of  surface  contamination  [Szuszczewicz  and  Holmes,  1975], 

density  fluctuations  [Holmes  and  Szuszczewicz,  1975;  Szuszczewicz  and  Holmes, 

1977]  and  magnetic  field  effects  [Szuszczewicz  and  Takacs,  1979].  Analysis 

of  approximately  25  characteristics  were  executed  over  the  F-layer  from  340- 

“3  -f 

560  km.  In  each  case  a  conversion  coefficient  a  [cm  ]/Ie(V  )  was  de¬ 
termined  so  that  the  I^fv"1")  profile  in  Figure  1  could  be  directly  scaled  to 

absolute  electron  densities.  This  procedure  yielded  a  =  (5.5  +  0.5)  10^ 

-3  -1 

electrons  cm  A  .  The  results  in  Figure  1  show  the  F-peak  at  375  km  with  a 

6  —  3 

maximum  density  of  1.3  (10  )  cm  (+  10%).  The  probe  measurements  also  re¬ 
vealed  a  number  of  maior  depletions  (6N  /N°  ^  90%)  distributed  throughout 

e  e 


the  F-region.  Each  of  the  depletions  has  been  shown  to  have  its  own  set  of 
smaller  scale  irregularities  [Szuszczewicz  et  al.,  1980;  Szuszczewicz  and 
Holmes,  1980]  with  the  most  intense  percent  fluctuations  occurring  on  the 
bottomside  gradient  near  270  km  (Region  "C"  in  Szuszczewicz  et  al.,  1980). 

The  pulsed  probe  data  provided  an  excellent  opportunity  for  comparison 
with  the  numerical  simulations  of  the  collisonal  Rayleigh-Taylor  (R-T)  insta¬ 
bility  at  intermediate  wavelengths  (e.g.,  Keskinen  et  al.,  1980).  Attention 
is  focused  on  the  depleted  structure  on  the  bottomside  F-layer  gradient  near 
270  km  which  is  believed  representative  of  the  mid-phase  development  of  the 
R-T  process.  (We  will  follow  Szuszczewicz  et  al.,  1980  and  refer  to  this 
domain  as  region  "C".)  An  expanded  view  of  this  domain  (Fig.  2)  shows  that 
it  is  not  a  single  macroscale  depletion  but  a  region  of  contiguous  large 
scale  structures  extending  over  a  total  altitude  domain  of  about  12  km. 
(Vehicle  velocity  in  region  "C"  was  2.4  km/ sec.) 

Typically,  computer  simulations  employ  several  values  for  the  zero- 
order  gradient  scale  length 


and  initialize  the  simulation  with  a  two-dimensional  perturbation.  In  the 
work  of  Keskinen  et  al.  [1980]  L  was  selected  to  be  5,  10  and  15  km  and  the 
perturbation  took  the  form([Chaturvedi  and  Ossakow,  1977]). 

6N  (x,y,t=0) 

— -  =  (10  H)  Sin  (k  y)  Cos  (k  x)  +  2(10  )  Sin  (2k  y) 

N  y  x  y 

e 
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* 

with  k  and  k  being  the  horizontal  and  vertical  wavenumbers,  respectively, 
x  y 

Both  k  and  k  were  set  equal  to  2ir/960  m  in  the  simulation.  In  addition, 
x  y 

the  computation  assumed  that  L  was  centered  at  300  km.  Under  actual  con¬ 
ditions  encountered  in  PLUMEX  I,  the  bottomside  F-layer  gradient  extended 
from  240  to  290  km. 

The  question  of  gradient  scale  length  can  be  studied  in  Figure  3  where 
it  is  shown  that  the  Dottomside  gradient  (encompassed  in  the 
105-125  sec  time  frame)  is  not  characterized  by  a  single  value  of  L.  Ex¬ 
perimentally,  L  was  computed  according  to  <1  >  (dl„/dy)-1  with  <1  >  and  the 

B  B  B 

altitude  derivative  determined  from  sliding  linear  fits  to  the  I  data  set. 

B 

In  region  "C"  (114s  ^  t  ^  122s)  L  is  seen  to  vary  between  2  and  10  km,  whereas 
adjacent  domains  (110s  £  t  <  113s  and  122s  £  t  <  126s)  can  be  characterized 
by  L  =  25  km. 

Previous  computer  simulations  [Keskinen  et  al.,  1980],  conducted  with 
L  =  5,  10  and  15  km,  gave  evidence  that  linearly  unstable  modes  saturate  by 
nonlinear  generation  of  linearly  damped  vertical  modes.  The  results  yielded 
one-dimensional  power  laws  (horizontal  and  vertical)  that  vary  with  a  spec¬ 
tral  index  (s  n  in  PN  <*  k  n)  between  2.0  and  2.5.  To  explore  this  result 

e 

within  the  context  of  region  "C",  and  to  provide  important  comparisons  with 

simulations  in  the  next  sections  ,  power  spectral  analyses  were  conducted 

over  sliding  intervals  of  2.4  km.  The  results,  presented  in  Figure  4,  show 

-2. 5 

that  the  dominant  behavior  is  k  '  over  the  range  k  =  2Tr/lkm  to  k  »  2tt / 2 5  m. 

—I  85 

The  k  "  behavior  at  t  =  116.001  sec  is  a  result  of  the  very  sharp  density 
gradient  (see  region  "C"  Fig.  2)  encompassed  by  the  domain  of  the  spectral 
analysis.  (We  note  that  the  spectral  analysis  is  executed  on  a  time  (fre- 
quencv)-domain  function  I  and  the  conversion  to  k-space  is  made  under  the 

D 

assumption  a  =  vf  \  where  v  is  the  vehicle  velocity.) 
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In  general  we  would  conclude  that  our  observed  large  depletions  and 
spatial  power  spectra  support  the  numerical  simulations  of  Kesklnen  et  al. 

(1980).  Further  testing  of  this  support  is  achieved  in  the  next  sections 
which  present  simulations  with  L  =  8  and  25  km  with  a  downward  drifting 
F-laver  model  that  is  more  in  keeping  with  the  actual  experimental  conditions. 
The  F-laver  time-history  can  be  important  since  unstable  modes  require  times 
in  excess  of  4,000  seconds  to  saturate... a  time  during  which  the  F-layer 
encountered  in  PLUMEX  I  drifted  downward  in  excess  of  40  km. 

MODEL  EQUATIONS 

The  two  fluid  equations  describing  the  Rayleigh-Taylor  instability  in  the 
presence  of  a  zeroth-order  horizontal  westward  electric  field  E  can  be  written, 
after  invoking  quasineutrality  (n^  »  n^  *=»  n) , 


dn 
_ 0 

at 


+  V  .  (n  V  )  =  -v  (n  -  n  ») 
aw  R  a  aO 


■(- 


V  x  B 

0  =  -T  Vn  -  en  I  -  Vcp  +  — - 10  +  E 

e  e  e  \  Y  c  -t 


) 


+  h  v)  h  '  +  “i  ( 


V.  x  B 

-Vcp  +  — — — +  E 
T  c  —o 


) 


+  mini5--miniVin^i 


(1) 

(2) 


(3) 


V  .  J  =  0  (4) 

J  =  n  (V.  -  V  )  (5) 

—  e  — i  -e 


where  a  denotes  electrons  (e)  or  ions  (i)  and  n,  cp,  V,  V  ,  V  ,  B  ,  and  g 

k  in  o 
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are  the  density,  electrostatic  potential,  velocity,  recombination  coefficient, 
ion-neutral  collision  frequency,  magnetic  field,  and  gravity  respectively. 

In  addition,  we  have  written  the  total  electric  field  E  =  -  V  cp  +  E^.  All 
other  symbols  retain  their  conventional  meaning.  Since  we  will  be  interested 
in  studying  intermediate  wavelengths  (X  ~  100m  -  1  km)  we  ignore  inertial 
and  pressure  effects  in  (1)  -  (5) .  Setting  Vcp  =  Vq>j_  ,+  m^g/e  [.Ossakow  et  al. , 
1979  J,  solving  algebraically  (2)  and  (3)  for  the  velocities  V  ,  _\/  ,  using  (4) 
and  the  ion  continuity  equation  (1),  we  find 


dn 

3t 


xz) 


V  n  =  -  v 
x 


R 


(n  -  n  ) 
o 


(6) 


V  . 
i 


CV.  nV.V  =  -  (gxz) 
in  i  i  c  — 


V  n  +  E  .  V  (v.  n) 
i  — o  ±  m 


(7) 


We  have  written  eq.  (6)  and  (7)  in  a  downward  moving  frame  whose  velocity  is 
2 

V  =  c  E  x  B/B  =*  10  m/sec  and  have  adopted  the  following  coordinate  system: 

L)  — O  — 

the  magnetic  field  B  is  in  the  z-direction  (north),  the  y-direction  (vertical) 

denotes  altitude  with  £.  =  -  gjy  and  the  x-axis  points  westward.  By  linearizing 

(6)  and  (7)  and  writing  6n  =  n  -  n  ,  cp,  ^  exp  U(k  x  +  k  y)  +  y,  tl,  k  .  B 

o  T1  x  y  k  —  — 

=  0 ,  and  |  Jk  |  L  »  1  we  f  ind 


Y 


k 


where  E  = 
— o 


R 


(8) 


From  (8)  we  note  that  we  recover  the  usual  collisional  Rayleigh-Taylor 
instability  if  E^  ->  0.  The  growth  rate  is  enhanced  (reduced)  by  an 
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eastward  (westward)  electric  field  E  .  We  will  be  concerned  only  with  a 

westward  electric  field  (E  =  E  x)  which  leads  to  a  downward  moving  F  layer. 

— c  — o  — 

NUMERICAL  SIMULATIONS 

Be  defining  n/  (x,v)  =  n  (x,v)/nQ(v),  (x,y)  =  .'x,y)/BL,  x  = 

x/L,  v'  =  y/L,  t'  =  v_t/L  where  n  (v)  =  Nq  (1  +  y/L)  (Nq  =  const.)  is  an 
equilibrium  solution  of  (6)  and  (7)  we  can  write  (6)  and  (7)  in  dimension-  y 
less  form  as  follows  (after  dropping  primes  for  clarity): 

v  dv-i  ->  So,  -  Sn  So. 

on  Sn  T1  cn  n  o  1  a  ,  ,  . 

St  ov  ox  Sx  Sv  n  S  Sx  H1  n 

o  y 


.2 

+  h+| 

fl  dn  +  1__ 

Sn  v 
_ °\ 

StPy 

[  1  dn 

Sc^! 

5x2 

,  2 

Sy 

n  Sy  n 

\  o 

Sv  I 

Sy 

n  Sx 

Sx 

A  Sn 
^2  n  Sx 

(10) 


where  p,  =  Lv„/c  and  B„  =  g/cv.  -  F,  /B  are  dimensionless  constants. 

Equation  (9)  was  integrated  forward  in  time  using  a  f lux-corrected  [Boris 

and  Book,  1973]  leapfrog-trapezoidal  scheme  [Zalesak,  1979]  while  eq.  (10) 

was  solved  for  the  self-consistent  potential  cp^  using  a  Chebychev  semi- 

iterative  multigrid  technique  [McDonald ,  1980 J.  The  computational  grid 

described  a  small  vertical  slice  oriented  in  the  east-west  direction  of 

the  bottomside  equatorial  F  layer  and  is  defined  by  64  x  64  points  with  a 

constant  mesh  spacing  of  Ax  =  Ay  =  1 3m  giving  an  altitude  and  east-west  extent 

of  960  meters.  Periodic  boundary  conditions  were  imposed  on  o/n  and  c,  in 

o  1 

both  the  x  (east-west)  and  y  (vertical)  directions. 

The  in  situ  rocket  data  of  Szuszczewlcz  and  Holmes  [1980]  do  not  reveal 
the  initial  conditions  from  which  the  observed  irregularities  developed. 
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Many  sets  of  initial  conditions  are  possible.  Tn  the  context  of  the  exr  ri- 

mental  observations  of  an  initially  upward  and  then  downward  moving  F  layer 

[Tsunoda,  1980a;  Szuszczewlcz  et  al.,  1980J,  we  make  the  reasonable  assumption 

that  the  irregularities,  whose  spatial  power  spectra  are  sampled  in  Region 

"C"  [Szuszczewicz  and  Holmes,  1980J  of  the  bottomside,  originate  at  a  higher 

altitude(s)  v  at  earlier  time(s)  and  are  convected  downward  bv  the  ambient 
'  o 

F  laver.  This  is  in  agreement  with  the  conclusions  of  Narcisi  and 
Szuszczewicz  (1980)  in  their  analysis  of  ion  composition  measurements  con¬ 
ducted  on  the  same  rocket.  From  previous  numerical  simulations  [Keskinen 
et  al . ,  1980]  of  the  intermediate  wavelength  collisional  Ravleigh-Tavlor 
instability  for  bottomside  gradient  scale  lengths  L  =  5  -  15  km  at  an  alti¬ 
tude  of  300  km,  we  found  that  a  well-developed  nonlinear  regime  could  be 

achieved  after  a  time  Zit  =*  10-15  v  ^  where  Y  is  the  maximum  linear  growth 

m  m 

rate.  We  define  the  altitude  of  instability  onset  y  to  be  the  position 

o 

of  the  bottomside  gradient  at  At  =■  10-15  y  ^  sec  prior  to  the  time  of 

rocket  penetration.  (The  rocket  penetrated  region  "C"  on  the  bottomside 

at  an  altitude  of  270  km.)  This,  of  course,  neglects  the  fact  that  y(y) 

will  be  decreasing  since  v  and  v  will  be  increasing  [see  eq.  (8)]  for 

a  downward  moving  F  layer.  In  other  wTords  we  wish  to  determine  an  altitude 

y^  -  15y  ^(y)Vp  =  270  km  where  v^  =  10  m/sec.  For  L  =  8  km  (25  km)  we  find 

y^  335  km  (365  km).  From  the  Jacchia  [1965]  model  neutral  atmosphere 

[Ossakow  et  al.  ,  1979],  (335  km)  =“  0.3  sec  \  v  (365  km)  =-  0.18  sec  \ 

V  (335  km)  =■  1  x  10  ^  sec  \  v  (365  km)  =“4.5  x  10  sec  ^  and  fi.  =  300 

R  R  i 

sec  ^ .  Coupling  these  considerations  with  the  probe  measurements  of  gradient 
scale  lengths,  two  different  numerical  simulations  were  made  using  L  =  8  and 
25  km.  Both  simulations  were  initialized  with  the  following  density  pro¬ 
file:  n(x,y,t  =0  )/N^  =  1  +  v/L  +  e(x,y)  where  e(x,y)  is  a  random  function 
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of  position  in  the  xy-plane  with  a  white  noise-like  spatial  power  spectrum  (no 

preferred  wavelength)  and  a  root-mean-square  amplitude  of  0.03. 

Since  the  real  space  dimensions  of  the  computational  grid  are  960  m  by 

960  m,  the  altitude  dependent  quantities  V1  (y)  and  VR(y)  do  not,  at  any  fixed 

time,  vary  appreciably  over  the  grid  since  their  scale  heights  are  much 

greater  than  the  grid  dimensions.  In  fact,  these  quantities  can  be  fit  very 

well  with  exponential  variation  in  y,  i.e.,  V  (y)  “  exp  (-y/L^)  and  VR(v) 

oc  exp  (-y/L  )  with  L  =“  55  km  and  L  =“  33  km  between  the  altitudes  of  250 
R  in  R 

and  550  km  using  a  Jacchia  j.  1965 J  model  neutral  atmosphere  [Ossakow  et  al. , 

1979].  1  ever,  over  several  thousands  of  seconds,  and  vR  will  change 

(increase)  in  a  frame  moving  downward  with  velocity  v^  =  1C  m/sec.  In  order 

to  compensate  for  this  effect,  we  convert  the  spatial  (altitude)  dependence 

of  v  and  vD  into  functions  of  time  by  making  the  substitution  y  -*■  -  v  t. 
in  R  d 

In  other  words,  during  the  course  of  the  simulations,  V^n(t)  =  (tQ)  exp 

(v  At/L  )  and  VD(t)  =  Vn(t  )  exp  (v  At/L)  where  At  is  elapsed  time  and 
L)  in  K  R  O  DR 

V.  (t  ),  v_(t  )  are  the  initial  values  depending  upon  initial  altitude.  In 
in  o  R  o 

this  simple  model  we  ignore  any  variation  of  the  bottomside  density  gradient 
scale  length  with  altitude.  We  now  present  the  important  nonlinear  aspects 
of  these  simulations. 

Figure  5  gives  an  isodensity  contour  plot  of  Sn(x,y)  =  bn/n^(y) ,  6n  = 
n(x,y)  -  nQ(y),  at  t  =  800  sec  for  L  =  8  km  and  shows  that  the  random  nature 
of  the  initial  perturbations  still  prevails.  The  contours  of  the  random 
initial  perturbations  in  §n(x,y)  at  t  =  0  sec  are  very  similar  to  Figure  5 
but  with  more  smaller  scale  structure.  Figure  6  displays  the  evolution  of 
6n(x,y)  at  t  =  2000  sec  where  some  vertical  elongation  and  steepening  can 
be  seen  together  with  small  scale  irregularities.  Figure  7  illustrates  the 
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perturbation  density  contours  at  t  =  4500  sec  where  further  steepening  and 
vertical  elongation  have  occurred.  This  late  time  density  configuration  is 
similar  to  recent  numerical  simulations  of  the  intermediate  wavelength 
collisional  Rayleigh- Taylor  instability  [Keskinen  et  al. ,  1980]  under  almost 
monochromatic  initial  conditions,  i.e.,  with  only  two  waves  initially  excited. 
The  maximum  percentage  depletion  (§n  <  o)  in  Fig.  7  was  56%.  Similar  density 
contour  development  and  late  time  percentage  depletion  were  found  using  L  = 

25  km  and  again  starting  from  random  initial  conditions.  These  large 
depletions  were  also  noted  in  the  in  situ  rocket  data  of  Szuszczewlcz  and 
Holmes  C 1980 3 • 

In  Figures  8a  and  8b  we  have  plotted  sample  one-dimensional  horizontal 

P(k  )  and  vertical  P(k  )  spatial  power  spectra  for  2ir/k  ,  2ir/k  between 
x  y  x  y 

approximately  100m  and  960m  in  the  nonlinear  regime  for  L  =  8  km  at  t  = 

/x*  2 

4000  sec.  These  power  spectra  are  obtained  by  integrating  |6n(k  ,k  )| 

x  y 

over  k^  and  k^,  respectively.  The  spectral  indices  obtained  from  these 
power  spectra  are  in  agreement  with  those  derived  from  bottomside  irregularity 
power  spectra  as  determined  from  rocket  observations  [Region  "C"  of 
Szuszczewicz  and  Holmes,  1980]  in  the  wavelength  domain  25m  -  1km  (see  Fig. 
4)  and  recent  simulations  ^Keskinen  et  al.,  1980].  Similar  power  laws  and 
spectral  indices  were  found  for  the  L  =  25  km  case  but  on  a  longer  time  scale. 

NONLINEAR  THEORY 

We  wish  to  interpret  the  evolution  of  long  wavelength  nlasma  irregu- 
larities  in  an  upward  or  downward  moving  equatorial  F  region  ionosphere  in 
terms  of  coherent  two-dimensional  mode-coupling  among  collisional  Rayleigh- 
Taylor  modes.  This  type  of  analysis  was  also  used  to  study  large  scale 
wavelengths  excited  by  the  gradient  drift  instability  in  the  equatorial 
electroiet  [Rognlien  and  Weinstock,  1974]  and  the  saturation  of  the  long 
wavelength  local  collisional  [Chaturvedi  and  Ossakow.  1977j  and  collisionless 
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L Hudson ,  1978]  Ravleigh-Tavlor  instability  In  equatorial  spread  F.  However, 

these  latter  papers  did  not  include  the  effects  of  a  convecting  (rising  or 

falling)  background  ionosphere  which  usually  accompanies  equatorial  spread  F. 

The  linear  growth  rate  vk  (eq.  8)  whicti  contains  the  altitude  dependent 

parameters  and  V^,  and,  by  inference,  the  nonlinear  growth  rate  will 

change  in  space  and  time  in  a  vertically  convecting  ionosphere.  Since  with 

vertical  drift  velocities  of  the  order  of  10  m/sec  the  F  region  ionosphere 

can  convect  over  a  distance  on  the  order  of  the  scale  height  of  u  and  u 

®  in  R 

on  time  scales  comparable  to  the  nonlinear  saturation  times  of  the  collisional 
Ravleigh-Taylor  instability,  this  effect  will  be  important. 

We  begin  our  nonlinear  analysis  with  eq .  (6)  and  (7)  which  can  be 
written 


«S  x  z  .  V  n  +  v  6n  =  £  V  6,  x  z  .  Vbn  (11) 

ot  o  j.  1  ~  o  R  B  x  l  — 

(c/B Q.)  V  .V.  n  V«5  =  (g  x  z/il.  +  cv.  F  /B)  .V6n  (12) 

ix  in  1  —  i  in— o 

where  we  have  made  the  separation  6n  =  n-n  .  We  note  that  the  addition  of 

o 

a  large  scale  electric  field  in  (12)  adds  only  a  linear  term.  The  con¬ 
vective  nonlinear  term  on  the  righthand  side  of  (11)  is  the  dominant  non¬ 
linearity  since  the  ratio  of  the  nonlinear  terms  in  (11)  and  (12)  is 


(c/Bii  )  (V  .V.  nV«h  )/(c/B)V  «S.xz  ,V6n 
i  x  m  1  x  1  — 


V ,  lit.  «  1 
in  l 
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for  two  dimensional  perturbations  {i^  and  6n.  This  allows  one  to  find  the 


potential  perturbation  &  for  arbitrary  density  perturbation  6n  usine  a 
linearized  version  of  (12),  i.e.. 


ecp 


ip 


6n 


(13) 


where,  keeping  the  same  notation  as  Chaturvedi  and  Ossakow  Ll977j  , 

2  2  2 

3  =  (8/v  ) (k  fl./k  c  )  with  c  =  T/m..  Equation  (11)  can  then  be  rewritten 

inxis  s  x 


S6n 

St 


'  \ 6" 


+  f 


x  z 


.  V6n 


(14) 


with  given  by  (8).  If  we  take  an  arbitrary  two-dimensional  perturbation 
of  the  form 


6n 

n 

o 


1,1 


sin 


(k  x 
x 


'UJt) 


cos 


k  y 
y 


(15) 


where  the  first  subscript  on  A  denotes  the  vertical  mode  number  and  the 
second  subscript  signifies  the  horizontal  mode  number,  we  find  using  (13) 
the  associated  potential  perturbation 


ecp 

— i  cos  (k^x  -  u>t)  cos  k^y  (16) 


Then  the  nonlinear  self-interaction  xz  ,  V  6n  generates  a  term 

of  the  form  n  a/2  A  ^  sin  2k  y  =  A  sin  2k  y  with  a  s  k^k  g/k^v,  which 
o  1,1  y  2,0  y  xy  in 

is  linearly  damped  by  recombination  using  (14).  This  linearly  damped  mode 
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14 


r 


V.  and  VD.  In  deriving  (18)  and  (19)  we  have  assumed  that  y(t)  and  a(t) 
are  weaklv  time  dependent.  For  0  (stationary  F  region),  we  recover 

the  results  of  Chaturvedi  and  Ossakow  [1977]  who  treated  y  and  a  as  space 
and  time  independent.  Since  a  ~  l/v^n,  one  observes  that  the  nonlinear 
coupling  decreases  (increases)  for  a  downward  (upward)  drifting  F  layer 
since  v  increases  (decreases).  Assuming  that  a  quasi-steady  state  is 
achieved  at  t  =  tg  with  dA^  ^/3t  “  g/dt  we  obtain  from  (18)  and 

(19) 


A2,0  (ts>  =  V1>1(ts)/2a(ts) 


1/2  k  L 

y 


A 

1,1 


(t  ) 
s 


<2  lY2,0<ts)|A2,0/“<ts>> 


(20) 


(<k2/kx,(VR(ts>vi„<ts)/8l-ky»'i 


(21) 


Since  Ao  ~(t  )  is  independent  of  v.  and  V,,  we  note  from  (20)  and  (21)  that 
t ,  u  s  m  k 

the  two  dimensional  spatial  power  spectra  of  the  density  fluctuations  in  a 
rising  or  falling  F  region  ionosphere  will  have  the  oblique  linearly  un¬ 
stable  modes,  e.g.,  A^  ^  being  favored  at  low  altitudes  (increasing  vR 

and  V,  )  with  the  linearly  damped  waves,  e.g.,  A  dominating  at  higher 
in  <2  ,U 

-2 

altitudes  (decreasing  V  and  V.  )  with  spatial  power  spectra  oc  k  .  Indeed 

R  in  y 

for  k  k  we  can  write 
y  x 
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so  that  at  300  km  with  L  =  25  using  a  Jacchia  [1965]  model  neutral  atmosphere 
[Ossakow  et  al . ,  1979]  we  find  ^  q  ~  ^  ^  ~  1%  for  k^L  =  2nL/X^  =-  50 
while  at  425  km  g/  ^  ~  10.  A  more  detailed  examination  of  the  static 
and  dynamic  solutions  of  (18)  and  (19)  will  be  reserved  for  a  future  report. 
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SUMMARY 


We  have  performed  analytical  and  numerical  simulations  studies  of  the 
collisional  Rayleigh-Taylor  instability  in  local  unstable  regions  of  downward- 
moving  equatorial  F  layers.  For  ambient  bottomside  plasma  density  gradient 
scale  lengths  L  =  8  and  25  km,  we  have  demonstrated  that  large  percentage 
relative  depletions  can  develop  on  time  scales  of  several  thousands  of  seconds 
from  purely  random  initial  conditions.  In  addition  we  have  shown  that  the 
one-dimensional  spatial  power  spectra  of  these  irregularities  in  the  vertical 
and  east-west  directions  conform  to  power  laws  k  n,  n  =  2-2.5  for  2ir/k 
between  100m  and  1  km.  These  results  are  in  good  agreement  with  our  recent 
in  situ  rocket  observations  of  intermediate  wavelength  (k:  25m  -  1  km) 
irregularities  equatorial  spread-F,  presented  here  in  this  paper,  complement 
previous  local  numerical  simulations  [Keskinen  et  al. ,  1980]  ,  and  lend 
further  support  to  the  belief  that  the  collisional  Rayleigh-Taylor  insta¬ 
bility  is  responsible  for  large  scale  size  irregularities  in  equatorial 
spread-F. 

Future  work  will  include  global  simulations  of  upward  and  downward 
moving  equatorial  F  layers  together  with  the  inclusion  of  neutral  wind 


effects. 
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Fig.  1  —  Relative  and  absolute  ionospheric  electron  density  profile  determined  in  situ  during  the 

occurrence  of  equatorial  spread  F. 
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TIME  AFTER  LAUNCH  [SEC] 

Fig.  3  —  Gradient  scale  lengths  on  the  bottomside  ledge  of  the  F-region  layer  shown  in  Figure  1. 
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8  0  0  SEC 


Fig.  5  —  Isodensity  contours  of  6n(x,y)/n0(y)  for  L  =  8  km  at  t  =  800  sec.  Solid  contours  denote 
8n/n0  >  0;  dashed  contours  denote  5n/n0  >  0.  Contours  are  evenly  spaced  with  y  axis  vertical  and 
x  axis  horizontal.  Tick  marks  and  numbers  denote  grid  point  locations  and  size  of  grid  in  meters, 
respectively.  Maximum  enhancement  (6n/n0  >  0)  and  depletion  (6n/n0  >  0)  are  +  4.8%  and  -  3.8%. 
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OICY  attn  &.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
[THACA,  NT  14850 

OICY  ATtn  0.  T.  FARLEY  JR 


ELECTROSPACE  SYStEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

OICY  ATTN  H.  LOGS TON 

OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 

OICY  ATtN  j.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 
OICY  Attn  C.  W.  PREMIE 

FORD  AEROSPACE  fc  COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  ALT0,  CA  94303 

OICY  ATTN  J.  T.  MA^INGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  OIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  8LVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  h.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OICY  ATTN  F.  RE  I BERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STRE£T  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATtN  DASIAC 
OICY  ATtN  DON  CHANDLER 
OICY  ATtn  TOM  BARRETT 
OICY  ATtn  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  Attn  WILLIAM  MCNAMARA 
OICY  ATtN  b.  GAMBILL 
OICY  ATTN  MACK  STANTON 

GENERAL  ELECtRIC  TECH  SERVICES  CO.,  INC. 

MMES 

COURT  STREET 
SYRACUSE,  NY  13201 

OICY  ATTN  G.  millman 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  80X  6770 
SANTA  BARBARA,  CA  93111 

OICY  ArTN  JOHN  ISE  JR 
OICY  AT^n  JOEL  GARBARINO 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBAMCS,  AK  99701 

CALL  CLASS  AtTN:  SECURITY  OFFICER) 
OICY  ATTN  T.  N.  QAVlS  (UNCL  ONLY) 

OICY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

OICY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

EL6CTR0NICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  of  ELECTRICAL  ENG  I  BERING 
URBANA,  IL  61803 

OICY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
uRBANA,  IL  60680 

CALL  CORRES  attn  SECURITY  SUPERVISOR  FOR) 
oicr  attn  K.  YEH 
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institute  for  defense  analyses 

400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  U.  M.  AEIN 
OICY  ATTN  ERNEST  SAUER 
01CY  ATTN  NANS  WOLFHARD 
OICY  ATTN  JOEL  SENGSTON 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  NA  01 7  JO 

OICY  ATTN  DONALD  NANSEN 


INTL  TEL  S  TELEGRAFN  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 


JAYCOR 

1401  C AMINO  DEL  MAR 
DEL  MAR,  CA  92014 

OICY  ATTN  S.  R.  GOLDMAN 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  M>  20810 

OICY  ATTN  DOCUfCNT  LIBRARIAN 

oicy  attn  Thomas  potemra 

OICY  ATTN  JOHN  DASSOULAS 

LOCKHEED  MISSILES  i  SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  0.  R.  CHURCHILL 

LOCKHEED  MISSILES  ANO  SPACE  CO  INC 
5251  HANOVER  Street 
PALO  ALT0,  CA  94504 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARD  G.  JOTNSON  OEPT  52-12 
OICY  ATTN  w.  L.  IFHOF  DEPT  52-12 

OMAN  SCIENCES  CORP 
P.  0.  BOX  7465 
COLORAOO  SPRINGS,  CO  80955 
OICY  ATTN  T.  FEAGHER 

L I FJCAB I  t  CORP 
10455  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  attn  IRWIN  JACOBS 

M.l.T.  LINCOLN  LABORATORY 
».  0.  BOX  75 
^EXINGTON,  MA  02175 

OICY  attn  OAVIO  M.  TOWLE 
OICY  Arx  P.  WALDRON 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  0.  CLARK 

MARTIN  marietta  CORP 
ORLANDO  DIVISION 
P.  0.  BOX  5857 
ORLANOO,  FL  52805 

OICY  ATTN  R.  HEFFNER 

MCOOAHELL  DOUGLAS  CORPORATION 
5501  80LSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.  W.  halprin 
OICY  attn  technical  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 

735  STA’E  STREET 

SANTA  BARBARA,  CA  93101 


OICY 

ATTN  P.  FISCHER 

OICY 

ATTN  W.  F.  CREVIER 

OICY 

ATTN  STEVEN  L.  GUTSCHE 

OICY 

at-tv  d.  sappenfielo 

OICY 

attn  r.  bogusch 

OICY 

ATTN  R.  HENDRICK 

QiCY 

ATTN  RALPH  KILB 

OICY 

ATTN  DAVE  SOWLE 

OICY 

ATTN  F.  FAJEN 

OICY 

ATTN  M.  SCHEI6E 

OICY 

ATTN  CONRAD  L.  LONGMIRE 

OICY 

ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 

P.  0.  BOX 

208 

BEDFORD,  MA  01730 

OICY 

AT^N  JOHN  morganstern 

OICY 

attn  G.  HARDING 

OICY 

ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  3LVO 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  ATtn  E.  C.  FIELD  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
518  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADORESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02175 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.  J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  10567 
OAKLAJC,  CA.  94610 

ATTN!  A.  THOMSON 


R  S  0  ASSOCIATES 

P.  0.  BOX  9695 

marina  DEL  RET,  CA  90291 

01CT  AT-N  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT  JR 
01CT  ATTN  ROBERT  F,  LELEVIER 
01CT  ATTN  WILLIAM  J.  KAR2AS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONA1D 
OICY  ATTN  r.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  8EDROZIAN 

RIVERSIDE  RESEARCH  INSTlrUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NT  10023 

01CT  attn  VINCE  TRAPANI 
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SCIENCE  APPLICATIONS,  INC. 

P.  0.  SOX  2551 
LA  JOLLA,  CA  92058 

Q1CY  A*TN  LEWIS  M.  -INSON 
01CY  ATTN  5ANIEL  A.  HA7R.IN 
01CT  ATTN  J.  SACHS 
01CT  ATTN  E.  A.  STRAKER 
01CT  ArTN  CURTIS  A.  SMITH 
11CT  AtTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  SOS’ON  POS'  ROAD 
SU08URY,  MA  01776 

OltY  ATTN  BARBARA  AOAMS 


SCIENCE  APPLICATIONS.  INCORPORATED 

1710  Goodridge  Dr. 
McLean,  VA  22102 
At ta ;  J.  Coctcayne 


Lockheed  Missile  &  Space  Co.,  Inc. 
Huntsville  Research  &  Engr.  Ctr. 
4800  Bradford  Drive 
Huntsville,  Alabama  35807 
Attn:  Dale  H.  Davis 


SRI  INTERNATIONAL 
355  RAVENSwOOO  avenue 
MENLO  PARK,  CA  99025 


OICY 

AT^N 

DONALD  NEILSON 

OICY 

ATT* 

ALAN  BURNS 

OICY 

AT-rN 

C 

Smith 

OICY 

A^N 

L  . 

L.  CCG6 

01CY 

AT-rN 

DAVID  A.  JOHNSON 

OiCY 

AT-rN 

WALTER  G.  CHESNUT 

Q1CY 

AT 

Charles  l.  rino 

OICY 

at-yy 

«AL 

TER  JAYE 

OICY 

AT^N 

M. 

BARON 

OICY 

ATN 

RAY 

’  L.  LEADABRANO 

OICY 

ATTN 

G. 

CARPENTER 

OICY 

ATTN 

G. 

PRICE 

OICY 

ATTN 

J. 

PETERSON 

OICY 

ATTN 

R. 

MAKE,  JR. 

OICY 

AT^N 

V. 

GONZALES 

OICY 

ATTN 

0. 

MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  S  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

01CY  ATtn  R.  <.  PLEBUCM 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

V1SIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WESt  industrial  park 

BURLINGTON,  MA  01803 

01CY  attn  Charles  Humphrey 
01CY  attn  J .  W.  CARPENTER 
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